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A general analytic procedure is developed for the post-Newtonian limit of /(i?)-gravity with 
metric approach in the Jordan frame by using the harmonic gauge condition. In a pure perturbative 
framework and by using the Green function method a general scheme of solutions up to {v/c) 4 order 
is shown. Considering the Taylor expansion of a generic function / it is possible to parameterize 
the solutions by derivatives of /. At Newtonian order, (v/c) 2 , all more important topics about the 
Gauss and Birkhoff theorem are discussed. The corrections to "standard" gravitational potential (tt- 
component of metric tensor) generated by an extended uniform mass ball-like source are calculated 
up to (v/c) 4 order. The corrections, Yukawa and oscillating-like, are found inside and outside the 
mass distribution. At last when the limit / — > R is considered the /(_R)-gravity converges in General 
Relativity at level of Lagrangian, field equations and their solutions. 
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I. INTRODUCTION 



The study of possible modifications of Einstein's theory of gravity has a long history which reaches back to the early 
1920s [ll-ll]- Corrections to the gravitational Lagrangian, leading to higher-order field equations, were already studied 



by several authors [4|-|6| shortly after General Relativity (GR) was proposed. Developments in the 1960s and 1970s 
(7H1 1| . partly motivated by the quantization schemes proposed at that time, made clear that theories containing only 
a R 2 term in the Lagrangian were not viable with respect to their weak field behavior. Buchdahl, in 1962 Q, rejected 
pure R 2 theories because of the non-existence of asymptotically flat solutions. 

In recent years, the effort to give a physical explanation to the today observed cosmic acceleration [l2Tjl4j has at- 
tracted a good amount of interest in /(i?)-gravity, considered as a viable mechanism to explain the cosmic acceleration 
by extending the geometric sector of field equations without the introduction of dark matter and dark energy. There 
are several physical and mathematical motivations to enlarge GR by these theories. For comprehensive review, see 

Other issues as, for example, the observed Pioneer anomaly problem 18] can be framed into the same approach 
[l|| and then, apart the cosmological dynamics, a systematic analysis of such theories urges at short scale and in the 
low energy limit. 

While it is very natural to extend Einstein's gravity to theories with additional geometric degrees of freedom, recent 
attempts focused on the old idea of modifying the gravitational Lagrangian in a purely metric framework, leading to 
higher-order field equations. Due to the increased complexity of the field equations in this framework, the main body 
of works dealt with some formally equivalent theories, in which a reduction of the order of the field equations was 
achieved by considering the metric and the connection as independent objects [20| . 

In addition, many authors exploited the formal relationship to scalar-tensor theories to make some statements 
about the weak field regime [2l| , which was already worked out for scalar-tensor theories [22j ■ Also a Post- Newtonian 
parameterization with metric approach in the Jordan Frame has been considered [23| ■ 

In this paper, we study the Post Newtonian limit of f(R) in the harmonic gauge. We are going to focus on the 
small velocity and weak field limit within the metric approach. In principle, any alternative or extended theory of 
gravity should allow to recover positive results of General Relativity. It will be very important to check at any level 
of our modified theory we can cover the outcomes of GR. 

The plan of the paper is the following. In the Sec]!]] we report the complete scheme of Newtonian and post- 
Newtonian limit of field equations for /(i?)-gravity and their formal solutions in the harmonic gauge condition. General 
comments about the mathematical properties of equations, their relative solutions (Gauss and Birkhoff theorem) and 
Minkowskian behavior of metric tensor are reported. In the Sec. IIIII we show the complete solutions (Newtonian 
and post-Newtonian level) when an uniform mass ball-like source is considered. The point-like source limit of the 
newtonian solution is considered and the compatibility of /(i?)-gravity with respect to GR is shown. Concluding 
remarks are drawn in Sec. IIVI 
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II. THE FIELD EQUATIONS UP TO POST-NEWTONIAN LEVEL 

Let us start with a general class of higher order theories given by the action 

A = J d 4 x^g~[f(R) + X£ m ] (1) 

where / is an unspecified function of curvature invariant R. The term C m is the minimally coupled ordinary matter 
contribution. In the metric approach, the field equations are obtained by varying ([1]) with respect to g^ v . We get 

/ / / 

H^v = f Rpu — -^Qtiv - f-,fj,u + g^Of = X T^ v (2) 

H = g afs H a(j = 2nf + f'R-2f = XT (3) 

Here, T p „ = ^^=^g ^ff is the the energy-momentum tensor of matter, while T = T° ' a is the trace, /' = , 
□ = ;0 . ;cr and X = 87TG 1 . The conventions for Ricci's tensor is R^ = R a ^av while for the Rienman tensor is 
R a Piiv = Fp u „ + ■■■■ The affinities are the usual Christoffel's symbols of the metric: V^ g — ^g^ [g a a,p+ gp<j, a — 9a0,<r)- 

The adopted signature is (H ) (sec for the details [HI)- 

The paradigm of post- Newtonian limit is starting from a develop of metric tensor (and of all additional fields in the 
theory) with respect to dimensionless quantity v. A system of moving bodies radiates gravitational waves and thus 
loses energy. This loss appears only in the fifth approximation in v. In the first four approximations, the energy of the 
system remains constant. From this it follows that a system of gravitating bodies can be described by a Lagrangian 
correctly to terms of order v in the absence of an electromagnetic field, for which a Lagrangian exists in general only 
to terms of second order. We thus find the equations of motion of the system in the next approximation after the 
Newtonian. 

To solve our problem we must start with the determination, in this same approximation of the weak gravitational 
field, of the metric tensor g^ v (for details see [25| ) 



l+g ( t ?(t, X )+g[t\t,x) + ... g { i\t, 



x 



^ l 9$\t, X ) + ... -^+ fl g ) (t,x) + ...J (4) 

The set of coordinates 2 adopted is x 1 * — (t, x 1 , x 2 , x 3 ). The Ricci scalar becomes 

R ~ RW(t,x) + R^(t,x) + ... (5) 
The n-th derivative of Ricci function can be developed as 

f n (R) ~ f n {R {2) + i? (4) + ...)- /"(0) + / ,l+1 (0)i? (2) + /" +1 (0)i? (4) + ^f n+2 (0)R {2)2 + ... (6) 
From lowest order of field equations (J2]) we have 



/(0) = (7) 

which trivially follows from the above assumption Q that the space-time is asymptotically Minkowskian. This result 
suggests a first consideration. If the Lagrangian is developable around a vanishing value of the Ricci scalar the relation 



1 Here we use the convention c = 1. 

2 The greek index runs between and 3; the latin index between 1 and 3. 
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([7]) will imply that the cosmological constant contribution has to be zero whatever is the /(-R)-gravity theory. This 
result appears quite obvious but sometime it is not considered in literature (27j . 
The Eqs. ([2j and © at 0(2) - order (Newtonian level) become 



' H$ = /'(0)i? t (2) - ^i?( 2 ) - /"(0)Ai?( 2 ) = XT t 



(o) 
// 



H 



(2) _ 



f(0)R 



(2) 



Z&RW +f"(0)AR^ 



*«-r(o)4 2 i=o 



(8) 



4 i?(2) = -3/"(0)Ai?(2) _ /'(0)i?(2) = AfT(°) 
where A is the Laplacian in the flat space, while at 0(3) - order become 



= /'(0)4 3) - f"(0)R (2 \u = XT^ (9) 

The solution for the gravitational potential g[f* /2 has a Yukawa-like behavior ([III]) depending by a characteristic 
length on which it evolves. Besides the Birkhoff theorem at Newtonian level is modified: the solution can be only 
factorized with a function space-depending and an arbitrary function time depending ([25]). Still more the corrections 
to the gravitational potential and the gravito-magnetic effects © are depending on the only first two derivatives of 
/ in R = 0. So different theories from the third derivative admit the same newtonian solution. 

Remembering the expressions of Christoffel symbols and using the following approximation for the determinant of 
metric tensor In y/^g ~ h\dtt ~ 9mm] + • ■ • , at 0(4) - order we have, 



-no) 



-f"(0) 



g { XmR (2 \n + Ai?W + , 9t (2) Ai?(2) + ff W Jj« _ i Vg £l . VR (2) 



\VRW\ 2 + RWARW 



= XT i 



(2) 



#(4) = _3/"(o)Ai?( 4 ) - f(0)i?( 4 ' - 3/"'(0) 



|Vi?^| 2 + i? (2) Ai?( 2 ) 



-3/"(0) 



ff (2) _ (2) „(2) _ i v , (2) _ (2) \ . y ff (2) _ (2) ff (2) 



= xt& 



(10) 



where V is the gradient in the flat space. Note that the propagation of Ricci scalar i?( 4 ) has the same dynamics of 
previous one. The complete knowledge of correction at fourth order for the rt-component of Ricci tensor fix the third 
derivative of / in R = 0. Also at this level there is a degeneracy of /(i?)-theory: different theories for only the first 
three derivatives admit the same gravitational field without obviously radiation emission. 

To complete the perturbative scheme and later find the solutions it needs to calculate the Ricci tensor components 
in j8]), ©, (fit)]) . After some calculus (pjj. psj) one obtains 



r r {2) - 

-"■it — 


I (2) 
2"tt,mm 




i? (4) - 


V 4) + 

2-ytt,mm ' 


1 (2) (2) 1 (2) (2) 1 ( 
2 gmn,mg t t,n + 29 mn 9tt,mn + 2#r 


R {3) - 


I„(3) _ 
2"ti,mm 


I (2) _ I (3) , 1 (2) 


d(2) 


±« (2) - 

2 -ytj^mm 


I (2) _1„(2) _1 (2) ,1 



1„(2) „(2) 



1 J2) 



,(3) 



(11) 



fin 



which represent the most general expressions without assuming any gauge condition. 
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A. The Newtonian limit of /(ii)-gravity 

We want to rewrite and generalize the outcome of ([25]) by introducing the Green function method (we remember 
that the Newtonian limit corresponds also to linearization of field equations) . Let us start from the trace equation. 
The solution for the Ricci scalar in the third line of (JSJ is 

i?( 2 )(t,x) = y£J d 3 x'g(x,x')T(°)(^x') (12) 

where m 2 = —jjtt^ and (?(x, x') is the Green function of field operator A — m 2 . 
The solution for g^\ from the first line of ([8} by considering that R^' — ^Ag^\ is 

fe^xj- 27r//(0) y^x |x _ x/| 4 Jrfx |x _ x/| 3m2 W (L6) 

We can check immediately that when / — > R wc find gf^(t,x) — > — 2G J d 3 x! ■ The expression (fl"3"|) is the 
"modified" gravitational potential (here we have a factor 2) for /(i?)-gravity. ^4 such solution which is the newtonian 
limit of f(R)- gravity is also gauge-free. 

Since we have a linearized version of field equations a such limit corresponds to one of Einstein equation and the 
linear superposition is satisfied. So the tt-component of energy-momentum tensor is, in this limit, the sum of mass 
energy volume density of sources: T t j = Ti a M a S(x — x a ) where <5(x) is the delta function. 

As it is evident the Gauss theorem is not valid since the force law is not oc |x| -2 . The equivalence between a 
spherically symmetric distribution and point-like distribution is not valid and how the matter is distributed in the 
space is very important (|29j). 

In this limit the geodesic equation is the Lagrangian of material point-like embedded in the gravitational field and 
its dynamics follows the Newtonian law: 

£ - ±V S <? (i ,x) (») 



B. The post-Newtonian limit of /(J?)-gravity in the harmonic gauge 

To simplify the expressions of the components of the Ricci tensor (|lip we can use the condition so-called of harmonic 
gauge: <7 po T^ CT = (see (25[). The Ricci tensor assumes the following simpler form: 



i? (2) 


= ±A g (2) 








p(3) 


= |A^f 








R {4) 


= |Aa£> h 


1 (2) (2) 
2~9mn9tt,mn 


I ( 2 ) _ 
2i>tt,tt 


il v^ 2) l 2 



(15) 



From the field equation ^ we find the general solution for g^ 



(3) 



5i 3) (*,x) = - 



X 



2vr/'(0) 



1 d 
6irm 2 dt 



(16) 



The choice of harmonic gauge enable us to solve with facility the equation © but we lose potential information 

(2) 

about the temporal dynamics of g\ t (t, x). A such knowledge is very important to obtain at least in perturbative 
approach some information about the Birkhoff theorem. By hypothesizing a perturbative approach (newtonian-like) 
we relegated inevitably eventual temporal dynamics only on the temporal variation of matter source. In fact in a such 
hypothesis of work the motion of bodies embedded in gravitational fields develops very slow with respect to motion 
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of matter. Then we have ever an instantaneous readjustment of spacetime. In other words the motion od bodies is 
adiabatic and it enables us to factorize the solution and with a time transformation we get a static solution. 



("2") 

From second line of ((U the solution for g\- follows 



1 

4tt 



, 3 ,i^(*,x') 



cfx' 



-^i?< 2 )(i,x) _ * J- / d 3 x'ii( 2 )(t,x') 
3m 2 67rm 2 |x| 3 J n 



9 I ,3 / d 3 x'i?( 2 )(i,x') 
27rm 2 |x| 3 7 n 



3m 



:i? (2) (t,x) 



.7. ; i 



(17) 



where fi| x | represents the integration volume with radius |x| (for the details see [3(|). By the solutions (THfl) . (|16l) . (fl7|) 
we can afhrm that it is possible to have solution non-Ricci-flat in vacuum: Higher Order Gravity mimics a matter 
source. It is evident from (|13|) the Ricci scalar is a "matter source" which can curve the spacetime also in absence 
of ordinary matter. Then it is clear also that the knowledge of behavior of Ricci scalar inside mass distribution is 
fundamental to obtain the behavior of metric tensor outside the mass. 

From the fourth order of field equation, we note also the Ricci scalar (R^ ) propagates with the same m (the second 
line of (ITUl) ) and the solution at second order originates a supplementary matter source in r.h.s. of ©■ The solution 
is 



i?(%x) = / dVe(x,x') j^T( 2 )(i,x') - 9 ix m (ty)R%Hty) - g^ n (ty)R%(ty) 



|Vtf.R< a > (t, x')| 2 + R (2) (t, x') A x ,i?< 2 ) (t, x') 



4 2) (t,x')-£L(*,xO 



•V x ,i?( 2 )(i,x') 



(18) 



where /i 4 = — gjwTgj ■ Also in this case we can have a non- vanishing curvature in absence of matter. The solution for 
glt \ from the first line of (|10p. is 



<4%x) = 



, l f XT$\t,x>) , i 



| Vi? (2) (t,x')| 2 + i? (2) (i, x') Ai?^ (t, x') 



1 

47T 



1 



67rm 2 



~ |X-X'| [ 2tt/'(0) 6tT^ 4 

^(t,x') 5 gL(t,x') - ff4 ( t 2) tf (t,x') - |V x , 5( (2) (^x')| 2 - i? (4) (t,x') - gW(t,x!)RW(t,x!) 
R^ _ ^(t,^(^ + „ ^ + ^ 



+ 5t (2) (t,x')A J R( 2 )(t ! x') +.g( 2 )(<,x') J R (2) , m „(t,x') - iv 5 ( n l(i,x') • Vi?< 2 )(t,x') 



(19) 



We conclude this paragraph by having shown the more general solution of field equations of /(-R)-gravity in the 
Newtonian and post-Newtonian limit assuming a coordinates transformation for the whose the gauge harmonic con- 
dition is verified. In the next paragraph we shall apply a such scheme to obtain the explicit form of metric tensor for 
a static and spherically symmetric matter source. 

III. THE SPACETIME GENERATED BY AN EXTENDED UNIFORM MASS BALL-LIKE SOURCE 

Let us consider a ball-like source with mass M and radius £. The energy-momentum tensor is (we are not interesting 
to the internal structure) 



T = p(x) 



(20) 
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where p(x) is the mass density and satisfies the condition g u u t 2 — 1, Ui — 0. Since ((4]) the expression ([20]) 
becomes 



T tt (t,x) ~ p(x)+p(x). 9t ( t 2) (<,x) = T<°^,x)+T t ( t 2) (t,x) 

t = P (x) = r(°)(*,x) 



(21) 



The possible choices of Green function, for spherically symmetric systems (i.e. Q(x,x l ) = £(|x — x'|)), are the 
following 



1 e" 



S(x,x') = 



47T |x — X'| 
i|*-x'| 



if 
if 



|x— x'| z |x— x'| 

with Ci + C*2 = — 4^. It notes, for any function of only modulus /i(|x|), that 



to 2 > 
m 2 < 



(22) 



7= /d 3 x'S(x,x>(x') = -^- / d|x'||x'| 2 M|x'|) / #'/ dfl' , = = = = (23) 

J 4ttJ Jo 7 ^Ixp + |x'| 2 - 2|x||x'| cosa 

where cosa = cos 9 cos 9' + siri0sin#'cos((^ — </>') and a is the angle between two vectors x, x'. In the spherically 
symmetric case we can choose 9 — without losing generality (the symmetry of system is independent by the angle) . 
By making the angular integration we get 



2m|x| 



d|x'||x'|/i(|x'|) 



e -m||x|-|x'|| _ e -m(|x|+|x'|) 



An analogous relation is useful also for the Green function of Newtonian mechanics |x — x' 



(24) 



, Mlx'l) _2tt 
Ixl 



= - ~|d|x'||x'| 



X — X' 



x - X — x - X 



Mlx'l) 



(25) 



A. Solutions at 0(2)- and 0(3)- order 

Supposing that m 2 > (i.e. sign[/'(0)] = — sign[/"(0)]) the Ricci scalar (fl2"|) . if we hypothesize a matter density 
asp(x) = H^-M), is 



1 - e-" lli (l + m£) 



sinh mlxl 



p -m\x\ 

e(e-|x|)-r 3 m 2 ^)-r^e(|x|-0 



(26) 



where 6 is the Heaviside function, F(x) = s mxcosil ^- smtimx and r g = 2GM is the Schwarzschild radius 3 . The 
solutions of (fll?j) . ([T^]) and (|T7|) . given the relations ([24]) and (j2"5)l , respectively are 



i 



2£ ' m 2 £ 3 



2£ 3 



" m «(l + m() sinhm|x| 
m 2 £ 3 m|x| 



e(^-|x|)-r 9 



1 

Ixl 



3 -m|x| 



9(|x|-0(27) 



9 t ( f (*,x) = 



(28) 



3 we have set for simplicity /'(0) = 1 (otherwise we have to renormalize the coupling constant X in the action JTJ). 
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<4 2) (i,x) 



(l + m£)e 



2£ 3m 2 £ 3 2f 3 



3m 2 £ 3 



2F(x) 



sinhmlxl 



e«-|x| 



x| 3m 2 |x| 3 3 \|x| m|x| 2 m 2 |x| 3 
2(1 + r7i£)e~" 1 ^ / sinhmjx 



-m|x| 



e(|x|-o U 



m 2£3 



\ m|x| 



f( x ; 



_2 2F(Q / 1 



6(£-|x|) 
3 



3 \|x| m|x| 2 m 2 |x| 3 



— m|x| 



e(|x|-o 



(29) 



(2) (2) 

For fixed values of the distance |x|, the solutions g\ t ' and g^' depend on the value of the radius £, then the Gauss 
theorem does not work also if the Bianchi identities hold [ll[. In other words, since the Green function does not scale 
as the inverse of distance but has also an exponential behavior, the Gauss theorem is not verified. We can affirm: the 
potential does not depend only on the total mass but also on the mass - distribution in the space. 

By introducing three metric potentials $(x), ^(x) and A(x) (the dimension is the inverse of length) we can rewrite 
(1271) and (HI) as follows 



' g%\t,y) =r g $(x) 



4 2) (t,x) = r 3 *( x )%+r 9 A(x)^ 
and with a fourth function, S(x), (the dimension is the cubic inverse of length) the Ricci scalar (|26|) is 



(30) 



^ 2 )(i,x) = r 3 S(x) (31) 

The spatial behavior ((21))) is shown in FIG.[TJ The metric potentials are shown in FIGs. [5J[3]and|U ^ * s interesting 
to note as the function $ assumes smaller value of its equivalent in GR, then in terms of gravitational attraction we 
have a potential well more deep. A such scheme can be intcrpretable or assuming a variation of the gravitational 
constant G or requiring that there is a central greater mass. These two affirmations are compatible on the one hand 
with the tensor-scalar theories (in the which we have a scaling of gravitational constant) and on the other hand with 
the theory of GR plus the hypothesis of the existence of the dark matter. In particular, if the mass distribution takes 
a bigger volume, the potential increases and vice versa. 



o.oo 

-0.02 
-0.04 
™ -0.06 
-0.08 
-0.10 
-0.12 



2 4 6 8 10 

M/f 

FIG. 1: Plot of dimensionless function £ 4 m -3 E for £ = ?n£ = .5 representing the spatial behavior of Ricci scalar at second 
order. In GR we would have H(x) = tt<9(£ — |x|). 




In the limit of point-like source, i.e. lim^ ~ l x l) — M5(x) and lim^o F(x) — 1, we get 
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FIG. 3: Plot of metric potential £m 1 \& vs distance from central mass with £ = m£ = .5. The dashed line is the GR behavior 
(similar to metric potential $). 



,2 e" 



,|x| 



' 1*1 

1 , le — l-l 



9 H T 3 |x| 



1 l 1 



9) |x| 3m 2 |x| a 3 I |x| m 2 |x| a l e 



(32) 



2 / 1 



2 |x|3 3 |x| T m|x| 2 ^ m 2 |x| 



3 lg-m|x| 



An important point has to be considered. The PPN-parameters 7 and j3, in the GR context, are intended to 
parameterize the deviations from the Newtonian behavior of the gravitational potentials. They are defined according 
to the standard Eddington metric (in the vacuum) 



r 9 



9tt — 1 |x| + 2 |x| 2 

k 9ij = Sij - iffiij 



(33) 
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2 4 6 8 10 



FIG. 4: Plot of metric potential £m X A vs distance from central mass with £ = m£ = .5. In GR a such behavior is missing. 

In particular, the PPN parameter 7 is related with the second order correction to the gravitational potential while (3 
is linked with the fourth order perturbation in v. Since the Gauss theorem is not verified for /(i?)-gravity, while the 
relations (J33J) satisfy it, we must consider the relations (I3"2"j) and not (j2"Tj) , ([23")) . After making the limit / — > R we can 
compare the results with (f3"3")l . Actually, if we consider this limit for ([3"2"| . we have 

' i?( 2 )(t,x) = 

< 9%\t,x) = -ft (34) 
k flgty.x) = -ft*« 

which suggest that the /(i?)-gravity is compatible with respect to GR. At last it is interesting to note that also in the 
case of extended spherically symmetric distribution of matter when we perform the limit / — > R the solutions (|27[) 
and (|29[) directly converge (in the vacuum) in solutions (|3~4"|) . demonstrating the validity of Gauss theorem in GR. 

Other important consideration is about the asymptotic behavior of /-theory with respect to GR. In fact increasing 
the distance from the central mass the gravitational field brings near to one of GR. A such rejoining is a normal 
consequence of the hypothesis (spherically symmetric system and request of Minkowskian limit). At last in FIG. 
Q] we report the spatial behavior of Ricci scalar (|26| approximating asymptotically the given value in GR. In fact 
hypothesizing a f(R)- theory the Ricci scalar acquires dynamics, and in the Newtonian limit, we find a characteristic 
scale length (m^ 1 ) on the which distance the scalar massive field evolves. Only for distances more great than m _1 
we recover the outcome of GR: R = 0. 

To conclude this section we show in FIG. [5] the comparison between gravitational forces induced in GR and in 
/(i?)-theory in the Newtonian limit. Obviously also about the force we obtained an intensity stronger than in GR. 

B. The oscillating Newtonian Limit of /(i?)-gravity 

If we consider m 2 < (i.e. sign[/'(0)] = sign[/"(0)]) from (|2"21 we can choose the "oscillating" Green function 



,. 1 cosmlx - x'l + sinmlx - x'l 
g(x,x) = -- : 35) 

4-7T X — X' 



The Ricci scalar (|12[) and the ti-component of at 0(2) order (|13|) become 



RW(t,x) 



6r„ 



1-H(0- 



0(£- |x|)- 2 r g m 2 G(0 



cos m x 



6(|x|-0 



(36) 



' 3 

.2! 



2£ 3 



2H{£) sinm|x| 



m 2^3 



1 2G(f) cosm|x| + sinm|x| 

IxT 3 W 



6(|x|-^037) 
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FIG. 5: Comparison between gravitational forces induced by GR and /(-R)-theory with £ = m£ — .5. The dashed line is the 
GR behavior. 

where G(£) = 3 mgcos ^ sinm ^ and#(£) = (l-m£) cos m£+(l+m£) sin m£, with the properties lim e _> G(0 = -1 
and lim£_>.o H (£) = 1. Since we have an oscillating Green function (it is not asymptotically zero) the "gravitational 
potentials" (|37[) at infinity are zero unless a constant value (lirria^oo 2r g mG(Q (sin ma — cos ma)). 

The spatial behavior of Ricci scalar (151)|) and metric component (13"T|) are shown in FIGs. |3]and[7] The considerations 



0.4 
0.3 

ra 

> 0.2 
0.1 
0.0 



10 20 30 40 

M/f 

FIG. 6: Plot of dimensionless function £ 4 m _3 E with £ = m£ = .5 representing the spatial behavior of Ricci scalar at second 
order in the oscillating case. 

of preceding subsection hold also for the solutions (1551) - (|37[) . The only difference is that now we have oscillating 
behaviors instead of exponential behaviors. The correction term to the Newtonian potential in the external solution 
can be interpreted as the Fourier transform of the matter density p(x). In fact, we have: 





lim / d 3 x'p(x')e" ikx = - limMG(|k|£) = M (38) 

Also in this case we conclude the section showing in FIG. [8] the comparison between gravitational forces induced in 
GR and in /(i?)-theory in the Newtonian limit. Obviously also in this last case we obtained a force stronger than in 
GR. 
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FIG. 7: Plot of metric potential C,m x $ vs distance from central mass with the choice ( = m£ = .5 in the oscillating case. 
The dashed line is the GR behavior. 




10 20 30 



FIG. 8: Comparison between gravitational forces induced by GR and /(_R)-theory with £ = m£ = .5 in the oscillating case. 
The dashed line is the GR behavior. 

C. Solutions at 0(4) - order 

The metric potentials and the function S(x), respectively denned in (|30)) and (f3l"j). satisfy the following properties 
with respect to derivative of coordinate Z-th in the matter 



4 we remember that |x| ; = |x| 1 x; 
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f S,,(x) = m2(1 + m?) e-™tF(x) Xl = S (x) 



J_ , (l+mg) 
F + 3^ 



l « F(x) 



1 (l+m£) - m f (m |x| 2 -f6) sinh mlxl+mlxl (m 2 |x| 2 — 6) cosh m|x| 

73 — e 



^jxT 5 



X; = <f (x)x; 



A/fx') = 2 ( 1 + m ^) e ~ m ^ ( 4m2 l x l 2 +9) si " h H x l-H x l(™ 2 l x l 2 +9)coshm|x| - = Anfxla;/ 

3 jIn (x) - 5 (x)ft n + 3(W) £ -mg ( m 2 |x| 2 + 3)sin h: n|x|-3 m |x|cosh m |x| ^ = ^ ^ + ^ 
^ $ lIn (x) = $ (x)& n + il+g^e-^ (- 2 ! x | 2 +3) sinh Hg~3m|x| cosh m |x| ^ = ^ (x)<W + ^ (x) ^ n 



and in the vacuum 

' S,«(x) 
$i(x) 

A,i(x) 



m (m|x|4-l) 



1 m|x| + l F(g)e- 

W + l x l a 3 



xi = $ (x)a;; 



m 2 |x| 5 



m 2 |x| 5 



m 2 |x| 5 



m 2 |x| 5 



a;/ = * (x)x; 
.T; = A (x)a:z 



S,« n (x) - So(x)g, w - TO (m |x l t 3TO|x|+3) F(0e- m|x| ^^ = S (x)«y i „ + S 1 (x)x ; a ;r> 



$,i„(x) = $o(x)<5/„ - 



3 , >» 2 |x| 2 +3m|x|+3 F(Qe~ 
|x|B + 1x1= 3 



X;X„ = $ (x)(5;„ + $i(x)xiX r , 



(39) 



(40) 



Obviously when we consider the physics in the matter or in the vacuum we choose the "right" quantities Sq(x), Si(x), 
$o(x), $i(x), * (x), A (x). 

The expression of Ricci scalar at fourth order, (ITS)) , is 



2 /-co 



2m|x| 7 



d|x'||x'N e 



'Ilir'lJ p - m ll x |-| x 'll _ p -™(l x l + l x 'l) 



,4 r 



,/|2 



S(x')^ + ^fS (x') 2 



3A(xQ + $0 ^-^^ +A °^]xf' 



So(x') + m 2 *(x')S(x') + 2 1 (x')A(x')|xf (41) 



from the which we note two contributes. The first one still depends on the quadratic term (oc R 2 ) in the action (JTJ, 
while the second one is related to cubic term (oc R 3 ). By introducing two functions £r(x) and 2//(x), the (|4Tj) is 
rewritable as follows 



i? (4) (i,x) =r g 



-A x ) + — =//(x) 



(42) 



An analogous situation is found for the tt-component of metric tensor at fourth order. In fact the (fT9")l becomes 
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5t ( t 4) (t,x) = ^^d|x'||x'||||x| - |x'|| - |x| - |x'||p(x')<f>(x') 



Rf ^Hxl{||x|-|x'||-|x|-|x'|}{i 



3*(x')+A(x') )$ (x')+ I *(x')+A(x') )$!(x')|x" 2 



3m 2 x 



— ' d|x'||x'H||x|-|x'||-|x|-|x'|U^- + A x ,S / (x') + S 1 (x')A(x')|x'| 2 



So(x') 



6A(x')+ (5*o(x')+3A (x') |x'| 2 



H(x') 



^ jT V|M{n - |x'|| - |x| - ! E "< V) - ix ' - l,X " J 



2m 2 $(x') + *(x') -3$ (x')-$i|x 



3m 4 



'|2 



A x ,5 7J (x') 
3m 2 



and by introducing other new functions $/(x), <I>jj(x) we have 



(43) 



9tt (*>x) = r 3 



*/W + TT $ ^(x) 



(44) 



It is useful to note that we have generally four contributes to gffl in (|4"3")l. The first one is induced by the non- 
linearity of the metric tensor even in our static spherically symmetric case. The product p(x)$(x) is not zero only 
in the matter but contributes to determination of tt-component in any point of space. The second one holds account 
of the induced contribution, by solution of previous order, to the determination of the tt-component of the Ricci 
tensor at fourth order. These first two terms are present also in GR. While the second two ones are derived from the 
modification of theory. In fact the third contribution depends on the addition of the quadratic term (oc R 2 ) in the 
action and finally the fourth one from the addition of the cubic term (oc R 3 ). 

The choice of free parameter (i (which is linked to third derivative of f{R)) is a crucial point in both the ex- 
pressions (|42l) and (|4"4")l to obtain the right behavior. From mathematical interpretation of Newtonian limit one 
has I/'" (0)| < |/"(0)| and if // > (i.e. sign[/'(0)] = - sign[/"'(0)], otherwise p 4 is not a length) we have 
to 4 //i 4 = |/"'(0)|/3/"(0) 2 , so we find the constraint < to 4 //! 4 < 1. In FIG. © we report the spatial behavior of 
(|42[) in the matter and in the vacuum, (.3 m 4 //i 4 ^ .9), showing that far from source we obtain a spacetime with a 
vanishing scalar curvature. At Newtonian level the Ricci scalar R^ is negative defined (ffi>| while at Post-Newtonian 
limit is positive defined. 

In FIG. (|10p we report the time-time component of metric tensor, g^\ on the same interval of values of to 4 //i 4 , 
although the behavior is quite insensitive to changes induced by the contributions of the cubic term in the Lagrangian. 
Besides we can observe an important analogy with respect the results of GR. In both cases we have a potential barrier, 
but for /(i?)-gravity it is higher (as in the Newtonian limit we found a deeper potential well). 



D. Solutions from isotropic coordinates to standard coordinates 

The found solutions of metric are expressed in isotropic coordinates and often for spherically symmetric problems 
is conveniently rewritten in standard coordinates (the usual form in which we write the Schwarzschild solution) . Here 
the relativistic invariant of metric Q is 



ds 2 = 



1 + r 9 $(x) 



$/(x) + — $ 7 j(x) 



dt 2 - 



1 - r 9 *(x) 



|dx| 2 +r g A(x) 



(x • dx) 



(45) 
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FIG. 9: Plot of dimensionless function C, 4 m 4 (2/ + zB.ii), representing the Ricci scalar at fourth order, where z = m 4 /^, 4 
and £ = m£ = .5. The spatial behavior is shown for .3 ^ z ^ .9 (solid lines) while the dotted line corresponds to 
R - jpU-^-theory. 




FIG. 10: Plot of dimensionless function C^ra 2 ($>i + z&u) (solid lines) where z = m 4 /fi 4 and of function l/2jx| 2 (dashed 
line). For z = (dotted line) we have the behavior of R — g^-ii 2 -theory. The solid lines are obtained for .3 ^ z ^ .9 and for 
s »<i .5. 

From spherically symmetric form of (|45|) it is convenient to replace the position x with spherical polar coordinates 
r, 9, (j) defined as usual by 

x 1 = rsin0cos</>, x 2 = rs'mOsincj), x 3 — rcoscj). (46) 
The proper time interval (|45l) then becomes 



da' 



1 + r g $(r) + r g 2 $j(r) + -^-$/r(r) 



di 2 



1 - rJ <f(r) + A(r) 



dr 2 



1 - r s <&(r) 



r 2 dft (47) 



where cZSl = d9 2 + sin 2 #d0 2 is the solid angle. To get the metric in the standard form it needs to impose a radial 
coordinate transformation 



1 - r g ^(r) 



1 ~2 

r = r 



and we have a new set of coordinates r, 9, (j). The metric (f4"7T) becomes 

dt 2 - 1 - r g [ *(f) + A(f) 



(is 2 = 



4 

m ~ 



1 + r ff $(f) + r/ ( $j(f) + — $ z/ (r 

A* 



^ I df 2 - f 2 df2 



(48) 



(49) 
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the expression previously desired. 

The explicit expression of P§j) is not displayed because the equation (|35]l can not be solved algebraically. However, 
this technical problem is overcome with the help of numerical methods when we are interesting to test experimentally 
the theory. In a further work that is taking place, we are considering the main experiments tested for GR and we 
want retest them with respect to the results obtained. 

IV. CONCLUSIONS 

In this paper, we have generally reformulated the Newtonian limit of /(i?)-gravity by applying the Green function 
method. Moreover the post-Newtonian limit has been studied in the harmonic gauge condition and the relative spatial 
behaviors (Yukawa-like and oscillating) of gravitational potential in the matter and in the vacuum have been shown. 
The Taylor expansion of a generic / has been considered obtaining general solutions in term of the derivatives up 
to third degree when an uniform mass ball-like source is considered. All metric potentials, however, depend strictly 
on the coupling parameters appearing indirectly in the Lagrangian of the theory (the derivatives of / in R = are 
arbitrary constants). 

A detailed discussion has been developed for systems presenting spherical symmetry. In this case, the role of 
corrections to the Newtonian potential is clearly evident. This means that one of the effects introducing a generic 
function of scalar curvature is to select a characteristic scale length which could have physical interests. 

Furthermore, it has been shown that the Birkhoff theorem is not a general result for /(i?)-gravity. This is a 
fundamental difference between GR and fourth order gravity. While in GR a spherically symmetric solution is, in any 
case, stationary and static, here time-dependent evolution can be achieved depending on the order of perturbations. 

Hypothesizing a nonlinear Lagrangian we obtained a gravitational attraction stronger than in GR. The hypothesis 
of dark matter is needed in GR to have more gravitational attraction. Here without hypothesis of alternative matter 
and without modifying the gravitational constant we have qualitatively the same outcome. This occurrence could be 
particularly useful to solve the problem of missing matter in large astrophysical systems like galaxies and clusters of 
galaxies. In fact dark matter could be nothing else but the effects that GR, experimentally tested only up to Solar 
System scales, does not work at extragalactic scales and then it has to be corrected. 

Then it is worth pointing out that /(i?)-gravity seem good candidate to solve several shortcomings of Modern 
Astrophysics and Cosmology. Taking into account the results presented, it is clear that only GR presents directly the 
Newtonian potential in the weak field limit while corrections appear as soon as the theory is non-linear in the Ricci 
scalar. 

In forthcoming researches, there is the intention to confront such solutions with experimental data in order to see 
if large self-gravitating systems could be modelled by them and if the experimental test of GR in the Solar System 
are compatible with /(i?)-gravity. 
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